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Abstract. Galaxies at high redshift with a strong Ly-a emission line trace massive 
O I star formation in the absence of dust, and can therefore be regarded as a prime 

I signature of the first major starburst in galaxies. We report results of the Ly-Q 

^ search within the Calar Alto Deep Imaging Survey (CADIS). With imaging Fabry- 



X 



Perot interferometer CADIS can detect emission lines in three waveband windows 
free of night-sky emission lines at 700 nm, 820 nm, and 920 nm. The typical flux 
detection limit for Ly-a emission redshifted into these windows, Fh™ 3 x 
lO^'^^Wm"^, corresponds to (unobscured) star formation rates of ^ lOM0/yr at 
z — 6. Candidate Ly-a-emitting galaxies are selected from the total emission line 
sample, which contains more than 97% of objects at 2: < 1.2, by the absence of flux 
below the Lyman limit (B-band "dropouts" ) , and the non-detection of secondary 
emission lines in narrow band fllters. We have detected 5 bright Ly-a-emitting 
galaxy candidates at 2 ~ 4.8, and 11 candidates at 2: ~ 5.7. For two of four 
observed Ly-a candidates, one candidate at z ~ 4.8, and the other at z; ~ 5.7, the 
emission line detected with the Fabry-Perot has been verified spectroscopically 
at the VLT. When compared to Ly-a surveys ai z < 3.5 even the upper limits 
set by our list of candidates show that bright Ly-a galaxies are significantly 
rarer at z <; 5 than the assumption of a non-evolving population would predict. 
Therefore we conclude that the Ly-a bright phase of primeval star formation 
episodes reached its peak at redshifts 3 < z < 6. 
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1. Introduction 

Important observational information about the epoch of galaxy formation can only be 
gained from the existence, number counts and properties of galaxies at very high redshift. 
In the past couple of years substantial progress has been made in discovering galaxies 
at redshifts z > 4. Several dozens of Lyman break galaxies have been detected at 2; ~ 4 
(e.g., Steidel et al. 1999). In addition, a substantial number of galaxies at 2: > 5 have 
been detected and verified spectroscopically (e.g., Dey et al. 1998; Weyman et al. 1998; 
van Breughel et al. 1998; Hu et al. 1999; Hu et al. 2002; Ellis et al. 2001; Dawson et al. 
2002; Ajiki et al. 2002; and Rhoads et al. 2002); one of these (van Breughel et al. 1998) is 
a radio galaxy at 2: = 5.2. Quasars have been identified out to z = 6.3 (e.g., Pentericci et 
al. 2002; Fan et al. 2001). However, there is a lack of systematic surveys for high redshift 
primeval galaxies. Many of the galaxies at 2; > 5 have been detected serendipitously. 

Moreover, as it has been pointed out, e.g., by Steidel et al. (1996), the strong metal 
absorption lines in the spectra of Lyman break galaxies indicate that they have already 
formed at least one generation of stars which has chemically enriched the systems and lead 
to dust formation, making Ly-a emission weak or undetectable. Therefore, the Lyman 
break galaxies do not show the properties expected for primeval galaxies. In contrast 
to the Lyman break galaxies, all but one of the galaxies found at 2; > 5 exhibit a very 
strong Ly-a emission line, in agreement with population synthesis models for primeval 
galaxies of Chariot & Fall (1993), which predict that a young, dust free, star forming 
galaxy should show strong Ly-a emission with (intrinsic) equivalent widths in the range 
of 5-25 nm. 

A direct probe of galaxy formation is to determine the number counts and redshift 
evolution of those nascent galaxies. Several systematic surveys aimed at discovering a 
significant number of 2 > 4.5 galaxies have now been proposed and are being carried 
out, e.g., Cowie & Hu (1998), Rhoads et al. (2000) and our survey, the Calar Alto Deep 
Imaging Survey (CADIS, Meisenheimer et al. 1998, 2003). At 2 = 4.5, a universe with 
standard parameters Hq = 70 km s~^ Mpc~^, Qm = 0.3, = 0.7 is only 1.2 Gyr 
old, corresponding to a look-back time of > 90% of the age of the universe. Any object 
observed at this or an earlier epoch should mark the early stages of galaxy formation. 

Unlike the Lyman break technique (e.g., Steidel & Hamilton 1992, 1993), which selects 
galaxies by the presence of young stars and intervening HI absorption, CADIS selects 
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galaxies by the extreme equivalent widths of their Lyman-a-emission {W\^^^^ > 10 nm). 
It has been demonstrated by Kudritzki et al. (2000) that the high Ly-a equivalent width 
of these galaxies can only arise in regions nearly free of dust, and therefore the Ly-a 
emission line is a good tracer of star formation from primordial material. Unlike similar 
programmes at 8—10 m telescopes, CADIS aims to find the most luminous starbursts, 
with SFR > 10 Mo/year, which most likely led to the formation of the galactic bulge. 
Even in the age of IG-m telescopes, only these luminous starbursts are accessible for 
detailed spectroscopic follow-up. 

This paper is structured in the following way: In section 2 we describe how we select 
Ly-a galaxy candidates at z > 4.7 with CADIS. In section 3 we present the spectroscopic 
follow-up of some Ly-a candidates. Finally, in section 4 we discuss the comparison be- 
tween observed and theoretical number counts of Ly-a-emitting galaxies up to z w 6. 

2. Lyman-a candidates from the CADIS survey 

CADIS combines a moderately deep multi-band survey (10 a limit, Riim = 24) with 
a deep emission line survey employing an imaging Fabry-Perot-Interferometer (Fum = 
3 • 10~^°Wm~^), which covers three waveband windows essentially free from atmospheric 
OH emission lines: A (A = 702 ± 6nm), B (A = 819 ± 5nm), and C (A = 920 ± 8nm). The 
classification and analysis of emission line galaxies from CADIS are outlined in Hippelein 
et al. (2003), Meisenheimer et al. (2003). Therefore, in this section we will present only 
those aspects which are relevant to the search of Ly-a galaxies using the CADIS (emission 
line) survey. 



Fig. 1. Night-sky emission spectra. The windows A, B, and C which are searched by CADIS 
are essentially free from atmospheric OH emission lines. 

The optimal detection of a weak emission line galaxy superimposed on the bright 
night sky is reached when the instrumental resolution SX is adapted to the expected 
line width. Assuming Av » 300kms~^, a typical virial velocity of a Milky Way- like 
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galaxy, yields SX = (1 + z)(Ati/c)Ao = 0.9 nm for a Lyman-a emission line at 2; ^ 6. 
On the other hand, the detection probability increases with the observed volume AV ~ 
AzAfi ^ AAA51, where Ail denotes the solid angle surveyed. Using a Fabry-Perot (FP) 
Etalon with AA = 2nm, which is tmied in 9 steps and covers an interval of 12 — 15 nm, 
CADIS provides an optimal compromise of both sensitivity and search volume. The three 
wavelength intervals (A, B, and C, see Figure 1) which are scanned by CADIS correspond 
to Lyman-a rcdshifts of z = 4.8, 5.7, and 6.6, respectively. The field size is Ail « 100 
arcmin^. Model predictions (Thommes & Meisenheimer 1995) of the galaxy abundance 
down to a line flux limit Funo > 3.0 • 10^^" Wm^^ (which can be reached with the 
Calar Alto 2.2- and 3.5-m telescopes within a few hours of observing time) made clear, 
before the project had even started, that at most a handful of Lyman-a galaxies would 
be expected per field and wavelength interval in the case of an open universe, which is 
commonly assumed today (Hm — 0.3, f^A — 0.7). Accordingly we decided to search every 
wavelength interval in five fields. After six years of observations on Calar Alto, 90% of 
the planned Fabry-Perot data have been obtained. The Fabry-Perot wavelength intervals 
which have been almost completely analysed are given in Table 1. The analysis of interval 
C has been postponed for the moment because of technical problems (strong "fringing" 
in the CCD images), and low expectations of success on the basis of the present results 
at 2; < 6, which will be presented in this paper. 

Table 1. The analysed FP- windows in four CADIS fields. 



Field 


FP-window 


Center coordinates 


An 


Az 


name 


observed 


q(2000) 


5(2000) 


arcmin^ 


analysed 






(hms) 


C ' ") 






Olh 


B 


01 47 33.3 


02 19 55 


105 


0.095 


09h 


A 


09 13 47.5 


46 14 20 


98 


0.082 




B 








0.070 


16h 


B 


16 24 32.3 


55 44 32 


107 


0.095 


23h 


A 


23 15 46.9 


11 27 00 


103 


0.057 




B 








0.082 



The following results are based on the data in the FP-windows in four CADIS fields 
presented in Table 1. It should be noted that the analysis of 69% of window B, and only 
27% of window A has been finished (see Table 1), in terms of the final AO x Az: AO, final 
will be » SOOarcmin^, corresponding to five CADIS fields, Az final for each Fabry-Perot 
window will be Azfinai ~ 0.1. 

The current classification of emission line galaxies shows that more than 97% of 
the emission line galaxies are foreground galaxies {z < 1.2), in which we mainly detect 
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the hncs Ha, H/3, [O III] A 5007, and [Oil] A 3727. Therefore, the main challenge is to 
separate those foreground objects from the Ly-a candidates. We identify Lyman-a-galaxy 
candidates among the emission line galaxies by the following exclusion criteria: 

(1) There should be no flux below the Lyman limit of a Lyman-a galaxy, because of an 
intrinsic drop in the spectra of hot O and B stars at the Ly- limit (Chariot & Fall 1993), 
and absorption by the neutral interstellar medium either in the galaxy itself (Leitherer 
et al. 1995a, 1995b) or in Lyman-limit foreground systems close to the primeval galaxy. 
For z > 4.7 the Lyman limit lies at A > 520 mn. Therefore, no flux should be detected 
in the CADIS B filter {Xcent — 461 nm, FWIIM= 113 nm). Accordingly, we require that 
Fb < 2 ■ (Tb ("blue dropouts"), where Fb is the flux measured in the B filter, and as 
its true overall error which takes into account both calibration errors and the scatter of 
the flux between the individual images. From 614 emission line galaxies selected in four 
fields in the FP-windows given in Table 1 this criterion yields NNobiue = 85 galaxies, i.e. 
13.8% of the 614 emission line galaxies. 

(2) Every single remaining Ly-a candidate has to be carefully checked on the FP, 
pre-filter and R images because of the possible contamination by nearby bright objects. 
Candidates which are closer than about 3" to a bright object are most likely spurious, 
since the changing atmospheric conditions change the tails of bright objects affecting the 
photometry of close-by galaxies. Therefore, such spurious objects are removed from the 
list of candidates. After this step, 5.7% of the emission line galaxies remain. 

(3) Foreground galaxies, for which we detect one of the prominent emission lines 
Ha, H/3, [O HI] A 5007, or [OH] A 3727 in the Fabry-Perot, can be excluded as soon a 
secondary line is detected in one of the veto filters, leaving 4.7% of the galaxies. 

(4) The distinction between a (rare) galaxy with bright [O II] A 3727 emission but 
weak restframe UV contimium, i.e., un(letec;table in the CADIS B filter, and a Lyman- 
a-galaxy is ambiguous based on the CADIS veto filters, because no line bluewards of 
[O II] A 3727 can be detected. Instead, we have to use the entire spectrum for this decision. 
Lyman-a galaxies at high redshift should distinguish themselves by a continuum step 
across the Lyman-a line, due to absorption from neutral hydrogen in the Ly-a clouds 
and Ly-limit systems along the line of sight (Madau 1995). Therefore, the candidates 
which remain after step (3) are classified in two categories: likely Ly-a, if they show 
almost no continuum on the blue side of the Ly-a emission line, e.g., no significant fiux 
in the R filter, and likely [O II] A 3727 objects otherwise. However, we need spectroscopic 
follow up observations of candidates belonging to both categories, since it cannot be 
excluded for sure, that some Ly-a galaxies fall between the likely [O II] A 3727 objects. 
After this step, 16 Lyman-a-candidates (classified as most likely Ly-a) remain (see Table 
2), i.e., 2.6% of the emission line galaxies we found in four CADIS fields. In addition we 
are left with 13 likely [O II] A 3727 candidates. 
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Table 2. The current hst of 16 CADIS Lyman-a-candidatcs. 

" The line of the Ly-a candidate has been verified by spectroscopic follow-up. 

" The line of the Ly-a candidate has not been verified by spectroscopic follow-up. 

'' Probable quasar. 



FPl 


Field 


Nr 


z 


Fline(Wm-2) 


A 


09h 


21556 


4.734 


5.7-10-2° 




23h 


34751 


4.772 


5.0-10-2° 




23h 


34105 


4.793 


4.2-10-2° 




23h 


50707" 


4.801 


4.0-10-2° 




09h 


38114 


4.741 


3.7-10-2° 


B 


23h 


40663 


5.732 


6.4-10-2° 




23h 


23836 


5.705 


5.7-10-2° 




23h 


28548 


5.694 


5.6-10-2° 




zon 


4o (4o 


c '700 
0.106 


0.1-10 




23h 


9324 


5.730 


4.8-10-2° 




16h 


3171 


5.746 


4.3-10-2° 




Olh 


3238" 


5.732 


4.1-10-2° 




23h 


45065 


5.735 


4.1-10-2° 




16h 


2314« 


5.694 


3.4-10-2° 




Olh 


27927" 


5.677 


3.1-10-2° 




Olh 


28090" 


5.681 


2.9-10-2° 



3. Spectroscopic observations 

Due to the low abundance of Ly-a galaxies (see section 4), spectroscopic follow-up at 
large telescopes constitutes an important part of the survey: First, we have to search close 
to the detection limit, which makes contamination by noise nonnegilible, and second, even 
very rare and unlikely contaminants, like distant supernovae or other transient objects, 
reach surface densities comparable to that of the Ly-a candidates, and can be mistaken 
as emission line galaxies in our observations of the FPI scans which are spread over years 
in some fields. 

The first goal of the spectroscopic follow-up of likely Ly-a galaxies with 8-m class 
telescopes is therefore to verify the emission line. Furthermore, the line shape and the 
continuum blue- and red-wards of the emission line could allow us to decide between Ly-a 
and [O II] A 3727: If the line has been verified, but the resolution is not good enough to 
identify the emission line, a second step is the clear confirmation of the line using higher 
resolution grisms, in order to see if the line shows the asymetric profile expected for Ly-a, 
or the double [O II] AA3726,3729 line with separation AA - (1 z), where AA = 2.75A . If 
galaxies with bright Ly-a line at ^; > 4.7 are confirmed, we plan a detailed study of the 
corresponding galaxy as third step. 
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Thus, confirming Ly-a candidates by spectroscopy is a slow process. Nevertheless, the 
verification or non- verification of a line for every single Ly-a candidate is very important 
in order to set more robust upper limits to the number counts of galaxies at high redshift, 
which can set stringent constraints on theoretical models (sec section 4). 

Spectroscopic follow-up observations of Ly-a candidates in the Olh- and 23h-fields 
were obtained in the summer and autumn of 2001, using FORS 2 at the VLT. The 
slitmasks contained 1.0" wide slits of length between 10" and 20". The position angle 
of each mask was chosen such that the number of Ly-a candidates in the 6.8' x 6.8' 
FORS field of view is maximized. Regions of the slitmasks not devoted to primary targets 
were used to obtain spectra of low metallicity emission line galaxies. The results of the 
observations of these low metallicity emission line galaxies will be presented and discussed 
in another paper (Maier et al. 2003, in preparation). 

Depending on the distribution of candidates, we used two different grisms, the lower 
resolution 3001 grism, which gives a spectral resolution of about 1.2 nm, at 800 nm, for 
slitlets 1" wide, and the 600 RI grism, which gives a higher spectral resolution of about 
0.8 nm, at 800 nm, for slitlets 1" wide. Three z 5.7 Ly-a candidates, Olh-3238, Olh- 
28090, and Olh-27927, were observed for 215 min using the 3001 grism with F0RS2, 
and one z « 4.8 Ly-a candidate, 23h-50707, was observed for a total of 150 min using 
the 600 RI grism with F0RS2. Fluxes were calibrated using multiple observations of 
the spectrophotometric standard stars LTT 7379, EG 274 and LTT 7987 (Hamuy et al. 
1992,1994). Seeing varied between 0".5 and 1".2. 

For two of the four observed Lyman-a candidates we verified the emission line, for 
Olh-3238 at 819.0 ± 0.3 nm, and for 23h-50707 at 705.7 ± 0.3 nm (see Figure 2 and Figure 
3). No continuum and no additional emission lines are seen on the VLT spectra, and, 
according to the CADIS measurements, these two objects satisfy the criteria (l)-(4) 
for Ly-a candidates. Therefore, we conclude that the two emission lines are very likely 
Ly-a fines of high redshift galaxies, Olh-3238 at z = 5.735 ± 0.003 and 23h-50707 at 
z = 4.803 ± 0.003, respectively. 

The line of the other two observed Ly-a candidates, Olh-28090 and Olh-27927, is 
not seen on the VLT spectra. These are the Ly-a candidates with the faintest fluxes of 
our Ly-a candidates list. One explanation for the non-detection could be that the slit 
possibly missed (part of) the line emitting region of these two galaxies. This is indicated 
by some emission line objects on the same masks as the Ly-a candidates which show a 
spectroscopic determined flux smaller than the CADIS flux. Thus, the slit width of 1".0 
used for the VLT observations was possibly too narrow, and we may have thus measured 
only a fraction of the emission line flux of the Ly-a candidates. The reason for this is that 
our astrometry delivers rather accurate coordinates (better than 0".2) for objects with 
multiple detections (FP and continuum), but can be off by more than 0".5 for objects 
detected only in the Fabry-Perot (Ly-a candidates). On the other hand, the two Ly-a 
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candidates could be indeed spurious, since we expect a residual contamination of about 
50% in our list of Ly-a candidates. The reason for the non- verification of the lines has 
to be established by extending the statistics of spectroscopic follow-up observations. 

It should be noted that we could verify an emission line (at 820.9 ± 0.3 nm) also for 
the galaxy Olh-4616. This object, which shows no flux in the B filter and in the veto 
filters, passed the criteria (l)-(3), but was classified as a probable [O II] A 3727 galaxy at 
z ^ 1.2 according to criterion (4). No other lines are seen on the VLT spectrum, and a 
continuum on both sides of the line is detected. Therefore, the emission line is indeed 
likely to be [Oil] A 3727 at z = 1.202 ± 0.001. 

We calculate the absolute Ly-a luminosity of Olh-3238 and 23h-50707 from the ob- 
served fluxes measured in the FP scan (since the spectroscopic observations are subject to 
slit losses, see above). The Lya flux is /(Lya) = (4.1 ±0.8) x 10"^° W m'^ for Olh-3238, 
and /(Lya) = (4.0 ± 0.8) x IQ-^o W m'^ for 23h-50707. Using a Qq = 0.3, Qa = 0.7, 
and Ho = 70 km s~^ Mpc~^ cosmology, we obtain an absolute Lya luminosity of i(Lya) 
~ (14.5 ±3.0) X 10^5 W for Olh-3238, and i(Lya) ~ (9.5 ± 1.9) x lO^^ W for 23h-50707. 
Note that these Lya luminosities are the highest found among any 2^5 galaxies (lu- 
minosities are estimated by using the same cosmology): e.g., 3.3 x 10^^ W for the lensed 
galaxy HCM 6A at ^ = 6.56 (Hu et al. 2002); 6.1 x lO^^ W for SSA22-HCM1 at ^ = 5.74 
(Hu et al. 1999); 3.4 x lO^^ W for HDF 4-473.0 at z = 5.60 (Weymann et al. 1998). Taking 
the Ly-a to Ha ratio for Case B recombination and no dust (Brocklehurst 1971), to- 
gether with the Kennicutt (1983) conversion between Ha luminosity and star formation 
rate, the derived star formation rates for the two Ly-a galaxies are: (14.5 ± 3.O)M0yr~^ 
for Olh-3238, and (9.5±1.9)M0yr-i for 23h-50707. Kudritzki et al. (2000) found that the 
derived star formation rates depend very strongly on the assumption whether we witness 
a burst or continuous star formation, the SPR being higher in the case of a starburst. 
Therefore, the derived SFRs using Kennicutt (1983) empirical formula for the case of 
continuous star formation represent only lower limits for the SFRs. The true values of 
the SFRs in case of a starburst may be a factor of five higher. 

4. Comparison between observed and theoretical abundances of Ly-a 
emitting primeval galaxies 

The CADIS selection suppresses the contamination in the list of Ly-a candidates to less 
than 50% (see section 2). As a consequence, the list of CADIS Ly-a candidates allows us 
to set stringent upper limits for the density of Ly-a-emitting galaxies at z > 4.7. We can 
therefore discuss the evolution of the population of Ly-a emitting galaxies from z = 3 
to ^; = 6 by comparing the maximum abundance derived from the list of CADIS Ly-a 
candidates with the results of other systematic surveys at 3 < 2; < 6. 
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Figure 4 shows the cumulative number counts of Lyman-a-galaxies at z = 4.8 (left 
panel) and at z = 5.7 (right panel), presented as the total number of galaxies, N per deg^ 
per Az = 0.1, which are brighter than a certain observed flux, -F/i„i- In the cmnulative 
histogram of CADIS Ly-Q galaxy candidates, the upper edge represents the total number 
of CADIS candidates per deg^ and per Az = 0.1, brighter than the respective Fum, 
based on the list of Ly-a candidates from Table 2: e.g., at z = 5.7, in four fields with 
a total Afltotai = 413arcmin^ « 1/9 we find one candidate brighter than Fum = 
6.4 • 10"2°Wm-2, two candidates brighter than Fum = 5.7 • 10"2°Wm"2^ ^^^^ The filled 
squares additionally take into account the likely residual contamination by noise plus 
unidentified [O II] A 3727-emitting galaxies, and therefore represent our best-guess upper 
limit to the abundance of Ly-a galaxies. 

Additional abundance measurements, showed as filled triangles, are from Hu et al. 
(1998) at z = 4.5, and from Hu et al. (1999) at z = 5.7. Hu et al. detected two Ly-a 
galaxies at z = 4.52 in a 24 arcmin^ field, and one Ly-a galaxy at z = 5.74 using 
narrowband observations in a 30 arcmin^ field. Stars indicate preliminary results from 
the Large Area Lyman Alpha Survey (LALA) at z = 4.6 (Malhotra & Rhoads 2002) and 
at z = 5.7 (Rhoads et al. 2002). Malhotra & Rhoads (2002) found 157 Ly-a candidates 
at 4.37 < z < 4.57 in one field of 0.36 deg^. The arrow in Fig. 4 (left panel) indicates the 
upper limit of the number density of Ly-a-cmitters based on the list of these candidates. 
Rhoads et al. (2002) reported 18 Ly-a candidates at z w 5.7 in one field of ^ 0.2 deg^, of 
which three (out of four observed) candidates have been verified spectroscopically. Rhoads 
et al. infer a number density 39 < < 54 of Ly-a-emitters per deg^ per Az = 0.1 down 
to their detection limit (see Fig. 4, right panel). For comparison, we show in Fig. 4 the 
measured abundance of Lyman-a-emitting galaxies at z = 3.5 through emission- line- 
surveys from Hu et al. (1998, small circles), and Kudritzki et al. (2000, large circles; 
converted from z = 3.1). 

In order to parametrize the observed number counts at z = 3.5, we use the model 
by Thommes & Meisenhcimer (2003): It basically assumes (i) that Ly-a-galaxies mark 
the onset of star-formation in present-day spheroids, (ii) a rather short-lived Ly-a-bright 
phase due to rapid dust formation (^ lO^yrs), and (iii) a mass dependent formation 
history as given by the peak formalism. In order to fit the abundance at z = 3.5 we 
used model parameters which set the peak of Ly-a-emission for the bulge of the Milky 
Way around z ~ 6. The normalization at z = 3.5 is set by the abundance of present 
day spheroids together with a free parameter — the duration of the Ly-a-bright phase. 
For the present work, we decided not to fine tune the model in a way that a consistent 
description of the abundances at all redshifts 3.5 < z < 5.7 is obtained. This will be 
presented in a forthcoming paper. Here we rather would like to demonstrate a robust 
qualitative statement — namely that Ly-a-galaxies are less common a,t z = 5.7 than at 
z = 3.5. To this end we simply shift the model function for z = 3.5 to higher redshift 
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by taking into account both the higher liuninosity distance (shift to the right), and the 
smaller comoving volume AV/Az (shift to bottom). This "no-evolution model" is shown 
by the dotted line in Fig. 4. 

It is obvious from Fig. 4 that the CADIS upper limits both at z = 4.8 and z = 5.7 
fall short of such a non-evolution model, while at the faint end, Fune < 10~^° Wm~^, the 
present results could still be compatible with no evolution. 

Since our statistics at 2; ~ 4.8 is extremely limited, we will focus any further discussion 
on redshifts around z = 5.7. 

While the number counts of Hu et al. at = 5.7 (a single galaxy found by Hu et 
al. 1999) still seem compatible with the no-evolution model, the abundance of Ly-a- 
galaxies from LALA and the limits given by the CADIS candidates fall significantly 
short of the no-evolution prediction. The under-abundance of Ly-a emitting galaxies at 
z = 5.7 compared to ^ = 3.5 becomes even more obvious if one considers that the noise 
distribution of our survey predicts residual contamination of about 50%. Nevertheless, 
before drawing definite conclusions, one should consider whether selection effects could 
lead to an under-estimation of the number of Ly-a-emitting galaxies. Potential influences 
on the upper limits include: objects classified as probable [O II] A 3727 galaxies for which 
the emission line is actually Ly-a, the fainter surface brightnesses of high redshift objects, 
and the large-scale structure. The weakly confined bright end of the luminosity function 
at ^; = 3.5 may limit the comparison at the bright end. Therefore, we have to discuss 
how robust the CADIS upper limits are. 

The distinction between a galaxy with [O II] A 3727 emission, but no flux in the CADIS 
B filter, and a Lyman-a-galaxy is ambiguous (see step (4) in Section 2). Therefore, 
some of the probable [O II] A 3727 objects may turn out to be Ly-a galaxies, and some 
presumed Ly-a galaxies may turn out to be [O II] A 3727-emitting galaxies. However, any 
asymmetry in this mutual contamination has already been accounted for by assuming a 
conservative estimate of 50% contamination of the Ly-a candidates sample. 

Surface brightness depends on the apparent size of an object at redshift z, which does 
not change much between z = 3.5 and z = 5.7, and on the luminosity distance. Since the 
model luminosity function is converted to z = 5.7 taking the luminosity distance into 
consideration, most of the effect of the (1 -|- z)~'^ dimming towards higher redshift has 
already been accounted for. 

Large-scale structure can influence galaxy counts even at high redshifts (Steidel et 
al. 1998). Since we derive the number counts of Ly-a galaxies at z pa 5.7 by combining 
four CADIS fields (with a total area about 14 times the size of the field searched by Hu 
et al. 1998), large scale structure should average out. Nevertheless, it should be kept in 
mind that results from small fields, like those searched by Hu et al. (1998, 1999), could 
be affected by large scale structure. 
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Two points should be noted about the bright end of the himinosity function at z « 3.5. 
First, the observed Ly-a hne fluxes for the most huninous galaxies at z « 3.5 found by 
Hu et al. (1998) and Kudritzki et al. (2000) are on the order of lO^^^Wm^^ Taking into 
account the higher luminosity distance at z = 5.7 (compared to z = 3.5) these galaxies 
would have apparent line fluxes of ^ 3 x 10^^" Wm^^ at z = 5.7 (using a JIq = 0.3, 
ft\ = 0.7 and Hq = 70 km s^^ Mpc^^ cosmology); i.e., galaxies with intrinsically bright 
Ly-a emission, corresponding to observed line fluxes > 3 x 10^^°Wm~^ at z « 5.7 — as 
searched by CADIS — have not yet been foimd in the small fields searched at z « 3.5. 
Thus, the shape of the bright end of the luminosity hmction relies mainly on our model. 
Second, the number counts at the bright end of the luminosity function at z = 3.5, 
resulting from the brightest galaxy detected by Kudritzki et al. (2000, large circles), are 
higher than the results of Hu et al. (1998, small circles). Since large-scale structure can 
influence galaxy counts even at high redshifts (Steidel et al. 1998), the Kudritzki et al. 
number counts could be biased, because Kudritzki et al. searched only one field, whereas 
Hu et al. survey goes over two fields. Disregarding the brightest galaxy from Kudritzki 
et al. would result in a steeper luminosity function at the bright end at z = 3.5 and thus 
lessen the discrepancy at z = 5.7. 

In summary, we conclude that the luminosity function of Ly-a-bright galaxies declines 
between z ~ 3.5 and z = 5.7. This implies that the peak of the Ly-a-bright phase — 
i.e., the first formation of massive stars — is reached after z = 6. Although, in principle, 
both luminosity evolution or density evolution could account for this finding, we only 
can establish the decrease in density around 3 x 10^^"Wm^^, as the bright end of the 
Ly-a luminosity function is still not determined well enough at z = 3.5. 

The onset of massive star formation should contribute a substantial fraction of the UV 
photons which reionize the universe. Moreover, both bright starbursts (SFR > 10 M0/yr) 
and the formation of supermassive black holes indicate a strong concentration of baryons 
within the centers of the most prominent density peaks. Therefore, our result seems to 
confirm additionaly the finding from the analysis of the HI absorption in the spectra of 
SDSS quasars at z > 6 (Fan et al. 2002), that the reionization of the universe did not 
happen long before z = 6. 
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Fig. 2. Olh-3238, a likely Ly-a galaxy at « = 5.732. 

(a) Photometry in all 14 optical CADIS filters fitted by a continuum- model (left panel), and the 
Fabry- Perot measurements with an emission line fit to the observed flux data (right panel) . 

(b) The single emission line, verified with FORS2 at the VLT, using the 3001 grism. 

(c) The two dimensional VLT spectrum. 
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Fig. 3. 23h-50707, a probable Ly-a galaxy at z ^ 4.801. 

(a) Photometry in all 14 optical CADIS filters fitted by a continuum- model (left panel), and the 
Fabry- Perot measurements with an emission line fit to the observed flux data (right panel) . 

(b) The single emission line, verified with FORS2 at the VLT, using the 600RI grism. 

(c) The two dimensional VLT spectrum. 
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Fig. 4. Cumulative number counts of Lyman-a-galaxies, A*' deg~^ per Az = 0.1, which are 
brighter than a certain observed flux Fum, a,t z — 4.8 (left panel) and z — 5.7 (right panel). 
In the cumulative histogram the upper edge represents the number of candidates, while the 
squares take into account the likely residual contamination by noise. Additional values of the 
number density are from Hu et al. (1998) for z = 4.5, and from Hu et al. (1999) at z — 5.7. 
Stars show preliminary results from the Large Area Lyman Alpha Survey (LALA) at z — 4.6 
(Malhotra & Rhoads 2002), and a.t z — 5.7 (Rhoads et al. 2002). Open circles denote the 
measured abundance of Lyman-a-galaxies at z — 3.5 through emission-line-surveys from Hu et 
al. (1998, small circles), and Kudritzki et al. (2000, large circles). A model luminosity function 
from Thommes & Meisenheimer (2003), fitted to the values at z = 3.5, has been converted to 
the redshifts searched by CADIS assuming no evolution between z = 5.7 und z = 3.5 (dotted 
lines) . 



